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3  Introduction 


Commercially  available  display  devices  for  use  in  digital  mammography  systems  are  known  to  have 
less  image  quality  than  the  mammographic  transilluminated  film.  Performance  parameters  that 
have  been  recognized  as  poor  include  low  maximum  luminance,  reduced  luminance  range,  veiling 
glare  and  ambient  light  reflections,  insufficient  resolution,  and  phosphor  granularity  noise,  among 
others.  However,  it  will  bep  ossible  to  achieve,  with  an  electronic  display,  better  image  quality 
than  that  observed  with  radiographic  film.  Organic  polymer  light- emitting  displays  (OLEDs)  are 
all-solid-state  flat  panel  devices  that  emit  light  through  electron-hole  recombination  processes. 
Recently,  OLEDs  have  received  worldwide  attention,  because  of  their  high  luminous  efficiency  that 
can  be  achieved  at  room  temperature  with  low  driving  voltages.  A  display  device  that  delivers 
better  image  quality  than  the  “gold  standard”  for  mammography  digital  systems  would  allow  a 
much  more  rapid  and  decisive  deployment  of  fully  digital  solutions  for  the  detection  of  breast 
cancer.  The  use  of  a  display  that  could  achieve  a  maximum  luminance  higher  than  2,000  cd/m^ 
while  having  a  very  low  minimum  luminance  would  allow  the  radiologists  to  implement  image 
processing  techniques  that  benefit  from  the  extended  luminance  range.  Such  a  display  would  allow 
for  film-less  practice  without  sacrificing  diagnostic  performance. 


4  Body 

The  electronic  display  of  mammographic  images  remains  today  an  important  technological  barrier 
for  the  deployment  of  fully  digital  mammography  systems.  Extremely  good  spatial  resolution  in 
large  image  pixel  array  sizes  are  required  to  detect  and  characterize  small  lesions.  In  addition, 
good  contrast  resolution  is  needed  to  achieve  high  diagnostic  performance  when  subtle  details  are 
present  in  low-luminance  image  regions  [1,2].  During  this  project,  we  have  achieved  the  following 
objectives: 

1.  Improve  available  simulation  methods  for  OLED  modeling  by  incorporating  actual 
photon  luminescent  spectra  characteristics  and  optical  material  properties  of  OLED  materials 
used  in  the  University  of  Michigan  laboratory. 

2.  Optimization  of  OLED  optical  performance:  We  have  used  advanced  computational 
modeling  tools  to  predict  the  achievable  improvements  and  define  new  structures  by  studying 
the  effect  of  surface  and  material  properties  on  the  optical  performance  of  AM- OLEDs. 

3.  Simulation  of  the  OLED  bidirectional  reflection  distribution  function:  We  have 
modeled  the  reflection  signatures  of  OLEDs. 

4.  Simulation  of  the  OLED  luminance  emission  angular  distribution.:  We  measured 
and  modeled  the  angular  luminance  emission  distributions  of  typical  OLED  structures. 

4.1  Simulation  methods 

After  the  simple  analysis  of  OLED  efficiency  presented  in  Ref.  [3],  several  groups  have  investigated 
the  effect  of  light  transport  in  multi-layer  structures  [4-7].  The  method  we  developed,  is  based  on 
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a  Monte  Carlo  (MC)  approach  [8,9].  The  MC  method  makes  use  of  the  generation  of  photons  with 
random  direction  according  to  a  distribution  function  describing  the  nature  of  the  light  emission. 
In  this  analysis,  the  light  source  within  the  organic  polymer  layer  is  considered  isotropic  from  a 
single  point  situated  in  the  center  of  the  device.  The  method  is  described  in  detail  in  Ref.  [10]  (see 
Appendix). 


4.2  Optimization  of  the  device  efficiency 

We  measured  the  index  of  refraction,  absorption,  photoluminescence  (PL),  and  electro- luminescence 
(EL)  of  three  organic  polymers  (A,  B,  and  C)  with  peak  emission  in  a  different  region  of  the  visible 
spectrum.  The  OLED  structure  modeled  is  a  hetero-structure  OLED  described  by  [11]  with  an 
aluminum  cathode  electrode  and  a  transparent  anode  ITO  electrode  (160  nm,  refractive  index  1.8). 
The  organic  polymer  film  thickness  used  was  200  nm.  The  transparent  substrate  index  of  refraction 
was  1.5  and  its  thickness  900  fim. 
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Figure  1:  Extraction  efficiency  r]pe,  waveguiding  rj^a  and  absorption  rjab  for  different  transparent 
conductive  oxide  (TCO)  electrode  models  (see  Ref.  [12]  for  details). 


In  Fig.  1,  we  present  results  showing  the  change  in  the  values  of  T/pe?  Vabj  and  r]^a  with  the 
index  of  refraction  (thickness)  of  the  transparentcond  ucting  oxide  (TCO)  electrode  covering  the 
range  from  1.1  to  2.7  (10-200  nm).  We  find  that  both  the  thickness  and  refractive  index  of  the 
transparent  electrode  have  a  minor  effect  on  the  shift  of  the  SL  spectrum  with  respect  to  the  PL 
emission.  We  obtained  a  maximum  r]pe  for  TCO  thickness  and  refractive  index  of  about  170  nm 
and  1.85  respectively.  The  waveguide  modes  are  determined  by  the  geometry  of  the  OLED  stack 
that  defines  a  total  internal  reflection  angle  {octir)  with  respect  to  the  device  plane,  beyond  which 
all  photons  are  emitted  through  the  edges.  Since  our  source  of  photons  is  isotropic,  we  expect 
lower  Tjma  foi*  larger  devices  due  to:  (1)  increased  probability  of  reflection  and  scattering  going  into 
the  solid  angle  defined  by  ariR^  and  (2)  increased  absorption  in  the  organic  polymer  layer.  For 
all  organic  polymers  considered  in  this  work,  our  results  show  that  the  simulated  light  emission 
is  shifted  toward  the  longer  wavelengths  consistently  with  experimental  measurements  of  EL.  We 
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also  show  that  the  photon  extraction  efficiency  is  reduced  by  light  absorption  and  waveguiding. 
The  reduction  of  the  OLED  wave-guided  fraction  can  be  obtained  by  using  structured  substrates, 
or  by  tuning  the  refractive  index  of  the  layers  to  maximize  transmission  and  reduce  total  internal 
reflection  at  each  interface. 


4.3  Reflection  properties  of  OLEDs 


The  effective  physical  contrast  of  display  devices  is  affected  by  the  achievable  contrast  generated  by 
the  light-producing  or  light-modulating  elements,  and  by  reflections  of  ambient  light.  In  emissive 
displays  that  are  typically  highly  reflective  by  design,  the  control  of  ambient  light  reflection  is  critical 
to  achieve  high  contrast,  especially  when  large  luminance  dynamic  range  is  needed.  For  instance, 
in  high-performance  medical  imaging  cathode-ray  tubes  (CRTs),  the  glass  faceplate  is  darkened  to 
reduce  diffuse  reflections  and  is  coated  with  an  anti-reflection  film  to  reduce  the  specular  component. 
In  this  phase  of  the  project,  paper,  we  developed  a  computational  method  for  modeling  the  BRDF 
of  electronic  displays  that  relies  on  an  optical  Monte  Carlo  simulation  code  developed  to  study  light 
transport  processes  in  emissive  structures.  The  BRDF  calculations  are  based  on  a  Monte  Carlo 
(MC)  method  developed  for  studying  light  transport  effects  in  emissive  structures  [8].  The  BRDF, 
a  formalism  often  used  in  optics  [13],  is  defined  by  a  six- dimensional  function: 


BRDF{ei,(l)uOo,(l>o,X,p) 


dLo{0oy  p) 

dEi(9i^  (j)i^  A,  p) 


(1) 


where  A  is  the  photon  wavelength  and  p  is  the  polarization  of  the  incoming  light  beam.  The  BRDF 
has  units  of  sr~^.  The  angle  of  incidence  of  ambient  light  is  defined  by  while  {9oy(l>o)  are 

the  angles  that  define  the  direction  of  reflected  light. 

The  precise  evaluation  of  this  function  is  time  consuming  and  costly.  The  BRDF  can  be  mea¬ 
sured  with  a  goniometric  setup  with  fixed  light  source  and  variable  detector,  with  variable  source 
and  fixed  detector,  or  with  a  conoscopic  approach  where  the  directional  intensity  is  mapped  into  a 
two-dimensional  distribution  recorded  by  a  position-sensitive  planar  detector.  The  first  two  meth¬ 
ods  suffer  from  severe  dependence  of  source  and  detector  positioning  precision  which  has  to  be  less 
than  1°.  The  third  method  requires  expensive  instrumentation.  The  models  of  display  devices  are 
based  on  uniform  slabs  representing  multiple  material  layers. 

In  this  study,  we  consider  only  the  integrated  response  to  an  unpolarized  light  source  with  flat 
spectrum  in  the  range  from  400  to  800  nm.  Finally,  the  angular  distribution  function  contained  in 
equally  spaced  bins  from  the  Monte  Carlo  simulations  is  reduced  to  a  luminance  distribution  density 
in  the  range  of  [0,90°]  for  a  given  incidence  angle  with  a  factor  l/cos{ao)*  The  introduction  of 
this  factor  is  needed  to  directly  relate  the  distribution  density  in  angles  with  the  expected  angular 
variations  in  the  luminance  out  of  the  device  (L).  is  also  known  as  the  in- plane  or  two- 

dimensional  BRDF.  Eq.  2  shows  the  final  metric  computed  from  the  Monte  Carlo  simulations. 


dL{ao)/cos{ao) 

dE{ai) 


(2) 


Organic  light-emitting  displays  are  simulated  using  a  typical  structure  model  that  includes 
measured  optical  properties  for  the  organic  polymer  films  and  the  substrate  [10].  The  OLED  model 
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(a)  (b) 

Figure  2:  B^^{ao)  for  (a)  a  monochrome  CRT  with  geometrical  and  material  models  described 
in  [8],  and  (b)  for  an  OLED  device  with  geometrical  and  material  models  described  in  [10]. 


represents  a  hetero-structure  OLED  described  by  He  et  al  [11]  with  an  aluminum  cathode  electrode 
and  a  transparent  anode  ITO  electrode  (160  nm,  refractive  index  1.8).  The  organic  polymer  film 
thickness  used  was  200  nm.  The  substrate  index  of  refraction  was  1.5  and  its  thickness  900  ^m. 
Fig.  2  shows  for  OLEDs.  It  can  be  observed  that  the  specular  component  is  predominant 

{Rs  —  0.10,  Rd  =  0.06  sr“^).  The  results  reveal  that  the  specular  peak  has  a  significant  width. 


We  performed  preliminary  experimen¬ 
tal  measurements  of  the  OLED’s  BRDF 
using  a  highly  accurate  BRDF  setup  at  the 
National  Institute  of  Standards  and  Tech¬ 
nology  (NIST)  in  the  Flat  Panel  Display 
Measurements  Laboratory  (directed  by  Ed 
Kelley).  In  Fig.  3  we  show  a  comparison 
of  the  experimental  high  resolution  BRDF 
data  with  the  predictions  of  the  simula¬ 
tion.  In  summary,  we  have  showed  that  the 
BRDF  of  OLED  displays  based  on  poly¬ 
mer  films  backed  with  a  metallic  reflective 
electrode  is  highly  specular.  We  reported 
for  the  first  time  the  complete  BRDF  sig¬ 
nature  of  OLEDs,  showing  that  the  haze 
component  is  small,  with  a  characteristic 
width  of  2°.  This  intermediate  component 
of  reflection  complicates  the  measurement  of  the  specular  peak  since  luminance  meters  having 
slightly  different  apertures  will  produce  a  specular  reflection  component  with  large  errors.  The  ex¬ 
perimental  results  agree  with  the  simulation  predictions  well  within  the  experimental  uncertainties 
present  in  even  high  resolution  BRDF  measurement  setups. 


Figure  3:  Comparison  of  high  resolution  BRDF  mea¬ 
surements  with  simulated  results  for  typical  OLEDs. 
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4.4  Angular  luminance  emission  of  OLEDs 


We  investigated  the  light  emission  angular  distribution  of 
the  OLEDs.  We  found  that  the  angular  distribution  of  the 
luminance  is  pseudo-Lambert ian  (rather  than  Lambertian) 
with  increased  intensity  in  the  forward  direction,  Fig.  4.  It 
was  found  that  the  light  refraction  at  different  interfaces, 
and  the  mirror-like  nature  of  the  bottom  surface  lead  to  this 
forward-peaked  angular  distribution,  while  the  absorption 
represent  a  minor  effect.  Finally,  it  is  clear  from  our  sim¬ 
ulation  that  the  waveguiding  effect  cannot  be  neglected  in 
the  OLEDs  and  its  contribution  will  be  dependenton  sub¬ 
strate  thickness  and  refractive  index  values  (see  Appendix). 

5  Key  Research  Accomplishments 

In  this  project,  we  have  achieved  the  following  accomplishments: 

1.  Development  of  improved  simulation  methods  for  OLED  modeling  by  incorporating 
actual  photon  luminescent  spectra  characteristic  of  OLED  materials  used  in  the  University 
of  Michigan  laboratory. 

2.  Simulation  and  optimization  of  OLED  optical  performance  with  respect  to  color  shift 
of  photoluminescent  spectra,  reflection  properties,  and  angular  luminance  distribution. 

6  Reportable  Outcomes 

During  the  project,  the  simulation  code  was  improved  and  ported  to  Windows-based  computers  and 
a  variety  of  Unix  machines  (Sun  UltraSparc  10,  HP  Vectra,  etc).  We  generated  one  manuscript  that 
was  published  in  the  Journal  of  Applied  Physics  [10],  and  two  presentations  [?,14]  (see  Appendix). 
Another  manuscript  is  in  preparation  [15].  The  simulation  code  DETECT-II  is  available  to  other 
investigators  on  a  “as-is”  basis.  For  information,  contact  Dr.  Aldo  Badano  (agbOcdrh  .fda.gov). 


Figure  4:  Simulated  angular  lumi¬ 
nance  distributions  of  OLEDs  for 
different  emissive  polymers. 
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Proceedings  paper  presented  at  the  International  Display  Research  Conference,  Osaka,  Japan, 
2001. 

Reprint  from  Journal  of  Applied  Physics:  "MONTE  CARLO  ANALYSIS  OF  THE  SPECTRAL 
PHOTON  EMISSION  AND  EXTRACTION  EFFICIENCY  OF  ORGANIC  LIGHT-EMITTING 
DEVICES”. 
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Bidirectional  Reflectance  of  Organic 
Light-Emitting  Displays 

Aldo  Badano,  Shu- Jen  Lee,  Jerzy  Kanicki,  Edward  F.  Kelley  and  Robert  J.  Jennings 


Abstract — The  reflection  properties  of  a  display  device 
influence  the  available  contrast  and  affect  the  perception 
of  low-luminance  detail.  We  report  on  a  Monte  Carlo 
method  for  modeling  the  bidirectional  reflectance  of  multi¬ 
layer  structures  used  in  thin-film  organic  light-emitting  dis¬ 
plays.  The  results  show  a  predominant  specular  peak  along 
with  a  quasi-Lambertian  component  and  significant  haze. 

L  Introduction 

The  effective  physical  contrast  of  display  devices  is  af¬ 
fected  by  the  achievable  contrast  generated  by  the  light- 
producing  or  light-modulating  elements,  and  by  reflec¬ 
tions  of  ambient  light.  In  emissive  displays  that  are  typ¬ 
ically  highly  reflective  by  design,  the  control  of  ambient 
light  reflection  is  critical  to  achieve  high  contrast,  espe¬ 
cially  when  large  luminance  dynamic  range  is  needed.  For 
instance,  in  high-performance  medical  imaging  cathode- 
ray  tubes  (CRTs),  the  glass  faceplate  is  darkened  to  reduce 
diffuse  reflections  and  is  coated  with  an  anti-reflection 
film  to  reduce  the  specular  component.  Previously,  we 
studied  the  effect  of  display  reflectance  on  the  physical 
contrast  in  the  low  luminance  end  of  the  grayscale,  and  re¬ 
lated  human  visual  system  performance  to  the  maximum 
ambient  illumination  that  maintains  a  perceptually  linear 
luminance  scale  [1].  It  has  been  shown  that,  for  CRTs,  the 
display  reflectance  can  be  characterized  by  the  sum  of  two 
components:  specular  and  diffuse  reflections.  However,  in 
flat-panel  displays,  the  complete  description  of  reflectance 
is  more  complex  and  requires  the  evaluation  of  the  bidi¬ 
rectional  reflection  distribution  function  (BRDF)  [5, 1 1]. 

In  this  paper,  we  describe  a  computational  method  for 
modeling  the  BRDF  of  electronic  displays  that  relies  on  an 
optical  Monte  Carlo  simulation  code  developed  to  study 
light  transport  processes  in  emissive  structures.  We  com¬ 
pute  the  angular  emission  distribution  for  a  beam  incident 
into  devices  having  different  structures.  We  present  re- 
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suits  for  Lambertian  and  specular  surfaces,  for  CRT  mon¬ 
itors  and  for  organic  light-emitting  displays  (OLEDs). 


II.  Methods 

A.  Monte  Carlo  simulations 

The  BRDF  calculations  are  based  on  a  Monte  Carlo 
(MC)  method  developed  for  studying  light  transport  ef¬ 
fects  in  emissive  structures  [1].  The  MC  method  makes 
use  of  the  generation  of  photons  with  random  directions 
according  to  a  probability  distribution  function  that  de¬ 
scribes  the  nature  of  the  light  source  [3, 8].  In  this  analy¬ 
sis,  the  light  source  is  a  collimated  beam  or  a  collimated, 
distributed  planar  source,  incident  on  the  display  surface 
with  a  direction  determined  by  (0i,  (j)i)  with  a  given  wave¬ 
length  distribution  S'(A),  and  polarization  state  p.  The 
photon  histories  are  followed  through  a  sequence  of  events 
that  includes  absorption,  scattering  and  Fresnel  refraction. 
The  simulation  models  bulk  absorption  events,  optically 
thin  coatings  and  rough  surfaces.  The  outcome  of  each 
individual  event  is  dependent  on  the  photon  energy  and 
polarization,  and  on  the  material  optical  properties.  Bulk 
absorption  is  determined  by  sampling  the  probability  of 
photon  absorption  after  a  path  of  length  /  using  an  expo¬ 
nential  law  Pabs{J>)  —  1  -  where  PabsW  is  the 

wavelength-dependent  linear  absorption  coefficient.  At 
the  optical  boundaries,  an  analysis  is  performed  depend¬ 
ing  on  the  surface  type  and  material  properties  using  Fres- 
neFs  equations  with  polarization  dependance  [7].  When 
the  film  thickness  is  comparable  to  the  photon  wavelength, 
we  use  modified  Fresnel  coefficients  to  describe  the  inter¬ 
ference  effects  of  optically  thin  films.  The  reflection  and 
transmission  coefficients  are  inteipreted  as  probabilities, 
and  therefore  the  photon  is  either  transmitted  or  reflected. 
The  simulation  outcome  is  calculated  by  statistical  aver¬ 
age  of  the  fate  of  all  histories  according  to  the  desired 
quantity  to  be  evaluated  for  each  experiment.  The  angu¬ 
lar  distribution  of  all  photons  reflected  by  the  device  is 
constructed  in  bins  with  varying  resolution  to  obtain  sat¬ 
isfactory  statistics  per  bin. 
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B.  BRDF  model 

The  BRDF,  a  formalism  often  used  in  optics  [9],  is  de¬ 
fined  for  any  reflecting  object  as  the  ratio  of  differential  re¬ 
flected  luminance  dLo  to  the  differential  illuminance  dEi 
incident  on  the  surface.  In  this  work,  we  consider  the  re¬ 
flectance  from  the  display  to  be  shift-invariant,  or  inde¬ 
pendent  of  position  across  the  screen,  therefore  neglecting 
all  edge-related  phenomena.  The  complete  expression  is 
then  given  by  a  six-dimensional  function: 


BRDF{9i,(j)i,eo,<t)oA.v) 


dLoi^Oo^  (po-) 

dEi(^0iy  (f>if  A,p) 


(1) 


where  A  is  the  photon  wavelength  and  p  is  the  polariza¬ 
tion  of  the  incoming  light  beam.  The  BRDF  has  units  of 
The  angle  of  incidence  of  ambient  light  is  defined 
by  {6i,  (pi),  while  (9o,  (po)  are  the  angles  that  define  the  di¬ 
rection  of  reflected  light. 

The  precise  evaluation  of  this  function  is  time  consum¬ 
ing  and  costly.  The  BRDF  can  be  measured  with  a  go- 
niometric  setup  with  fixed  light  source  and  variable  de¬ 
tector,  with  variable  source  and  fixed  detector,  or  with 
a  conoscopic  approach  where  the  directional  intensity  is 
mapped  into  a  two-dimensional  distribution  recorded  by  a 
position-sensitive  planar  detector.  The  first  two  methods 
suffer  from  severe  dependence  of  source  and  detector  po¬ 
sitioning  precision  which  has  to  be  less  than  1°.  The  third 
method  requires  expensive  instrumentation. 

The  models  of  display  devices  presented  in  this  paper 
are  based  on  uniform  slabs  representing  multiple  material 
layers.  In  flat  panel  display  devices,  the  specific  details  of 
the  pixel  structure  influence  reflection  and  propagation  of 
light  within  the  different  layers.  In  this  model,  we  incor¬ 
porate  the  effect  of  device  features  at  the  pixel  level  into  a 
general  description  of  the  scattering  properties  at  the  sur¬ 
face  between  the  pixel  circuits  and  the  substrate.  With  this 
assumption,  the  light  scattering  processes  can  be  consid¬ 
ered  radially-symmetric,^  and  the  BRDF  can  be  described 
by  a  function  given  by 


B'{ai,aoA^  P) 


dLo(^^o)  A,  P) 
dEi{ai,  A,  P) 


(2) 


In  this  study,  we  consider  only  the  integrated  response 
to  an  unpolarized  light  source  with  flat  spectrum  in  the 
range  from  400  to  800  nm.  Therefore,  integrating  over  all 
A  and  p  for  a  given  light  source  described  by  S' (A)  and  p 
yields  a  simplified  function 


B^^iao)  =  dLo{ao)/dE^{ai)  (3) 

^This  assumption  is  not  true  for  emissive  displays  with  polarizer 
films  or  for  liquid  crystal  displays.  An  extension  of  this  work  to  in¬ 
clude  such  structures  is  in  progress. 
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Fig.  1 .  B{a)  for  a  surface  with  50%  specular  and  50%  diffuse  compo¬ 
nents.  The  specular  fraction  results  in  a  delta  function  at  the  specular 
angle  of  incidence  (in  this  case  1 0°).  The  correction  with  a  factor  equal 
to  cos  (a)  yields  a  horizontal  profile,  typical  of  Lambertian  surfaces. 


Finally,  the  angular  distribution  function  contained  in 
equally  spaced  bins  from  the  Monte  Carlo  simulations  is 
reduced  to  a  luminance  distribution  density  in  the  range 
of  [0, 90°  ]  for  a  given  incidence  angle  ai  with  a  factor 
l/cos(ao).  The  introduction  of  this  factor  is  needed  to 
directly  relate  the  distribution  density  in  angles  with  the 
expected  angular  variations  in  the  luminance  out  of  the 
device  (L).  P"‘(ao)  is  also  known  as  the  in-plane  or 
two-dimensional  BRDF.  Eq.  4  shows  the  final  metric  com¬ 
puted  from  the  Monte  Carlo  simulations. 

=  dLMcosM  (4) 

dJlf  \CXi  ) 

Let  us  consider  a  reflective  smface  for  which  half  of 
the  time  photons  are  reflected  specularly,  and  the  rest  of 
the  time,  undergo  diffuse  reflection  according  to  a  Lam¬ 
bertian  distribution.  For  this  case,  the  BRDF  is  known 
analytically  and  is  given  by  Eq.  5,  taking  into  account  the 
number  of  histories  (H)  and  the  number  of  angular  bins 
(/?), 

=  n,  m  “T" 

^(2  X  0.5  X  H)  ao  ^ 

The  factor  of  2  comes  from  the  summation  of  photons 
at  positive  and  negative  angles  with  respect  to  the  sur¬ 
face  normal.  Fig.  1  shows  the  output  angular  distribution 
B^^{ao)  of  an  ideal  Lambertian  surface  with  a  50%  spec¬ 
ular  component.  The  expected  flat  response  of  a  Lamber¬ 
tian  surface  is  seen  in  the  simulation  results.  For  this  case, 
the  number  of  angular  bins  provides  an  angular  resolution 
of  0.009°  (/?  --  10^),  therefore: 


Fig.  2.  Detected  photon  density  using  a  1 00  x  1 00  binning  array  cor¬ 
responding  to  a  high  resolution  monochrome  CRT,  simulated  using  a 
sensor  array  in  close  proximity  to  the  display  screen. 


Fig.  3.  J5^°(ao)  for  a  monochrome  CRT  with  geometrical  and  ma¬ 

terial  models  described  in  [1]  for  an  incidence  angle  of  10°  using 
H  =  10®. 


5  X  10®  ao  =  10°, 

10^  Oo  #  aj. 


(6) 


The  dimensionless  specular  reflection  coefficient  Rs  is 
approximately  given  by  the  ratio  of  the  raw  specular  peak 
amplitude  to  the  total  number  of  histories: 


Rs 


number  of  photons  in  the  specular  peak 
total  number  of  photon  histories 


(7) 


which  in  this  case  is  Rs  —  0.5.  Analogously,  the  diffuse 
reflection  coefficient  Rd  (in  nit/lux)  can  be  approximated 
by  the  ratio  of  the  fraction  of  non-specular  photons  over 
27r  to  compensate  for  the  emission  into  the  hemisphere: 


Rd 


1 

(3  27r 


=  0.02 


(8) 


The  corresponding  detected  photon  density  per  pixel  us¬ 
ing  a  lOOx  100  binning  array  at  the  exit  plane  is  shown  in 
Fig.  2.  The  in-plane  photon  density  is  peaked  at  the  center 
where  the  incident  beam  impinges  on  the  display. 


III.  Results 

We  present  results  for  incident  collimated  sources  at 
10°.  We  used  a  fixed  light  source  with  flat  spectral  dis¬ 
tribution  from  400  to  800  nm.  For  instance,  Fig.  1  shows 
the  output  angular  distribution  of  an  ideal  Lambertian  dif¬ 
fuser  surface  with  a  50%  specular  component.  It  can  be 
shown  that  the  specular  reflection  coefficient  Rs  is  ap¬ 
proximately  given  by  the  ratio  of  the  raw  specular  peak 
amplitude  to  the  total  number  of  histories,  and  that  Rd 
can  be  calculated  using  Eq.  8  [2]. 

To  model  typical  high-performance  medical  imaging 
CRT  monitors,  we  used  a  structure  validated  in  previous 
work  [3].  The  CRT  emissive  structures  consisted  of  an  ab¬ 
sorptive  glass  faceplate  of  13  mm  with  a  phosphor  layer 
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Fig.  4.  for  an  OLED  device  with  geometrical  and  material 

models  described  in  [4]  for  an  incidence  angle  of  10°  using  H  =  10®. 


having  a  diffuse  reflectance  of  0.90-0,99,  and  an  anti- 
reflective  coating.  Fig.  3  shows  the  specular  and  quasi- 
Lambertian  components  typical  of  CRTs.  The  value  of 
the  specular  coefficient  Rs  calculated  from  the  raw  spec¬ 
ular  peak  is  around  0.045,  and  Rd  is  around  0.02  sr~^. 
These  values  agree  well  with  previous  measurements  of 
similar  devices  [1].  The  decrease  in  luminance  with  angle 
can  be  associated  with  the  absorption  in  the  faceplate  and 
most  importantly,  with  the  non-Lambertian  nature  of  the 
reflections.  We  conclude  from  these  results  that  the  dif¬ 
fuse  reflections  (and  possibly  the  emissions)  from  CRTs 
are  more  forward-peaked  than  would  be  the  case  for  a  per¬ 
fect  Lambertian  diffuser.  These  results  are  consistent  with 
recent  observations  reported  by  Blume  [6]  where  the  lu¬ 
minance  of  CRTs  was  measured  at  different  viewing  di¬ 
rections  using  a  spot  photometer. 

Organic  light-emitting  displays  are  simulated  using  a 
typical  structure  model  that  includes  measured  optical 
properties  for  the  organic  polymer  films  and  the  sub¬ 
strate  [4].  The  OLED  model  represents  a  hetero-structure 
OLED  described  by  He  et  al  [10]  with  an  aluminum  cath¬ 
ode  electrode  and  a  transparent  anode  ITO  electrode  (160 
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Fig.  5,  Normalized  reflected  spot  profile  for  the  CRT  model. 
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Fig.  6.  Normalized  reflected  spot  profile  for  the  OLED  model. 

nm,  refractive  index  1.8).  The  organic  polymer  film  thick¬ 
ness  used  was  200  nm.  The  substrate  index  of  refraction 
was  1.5  and  its  thickness  900  /um.  Fig.  4  shows  B^^{ao) 
for  OLEDs.  It  can  be  observed  that  the  specular  compo¬ 
nent  is  predominant  (Rs  =  0.10,  Rd  ~  0.06  sr~^). 
The  results  reveal  that  the  specular  peak  has  a  significant 
width.  This  observation  is  consistent  with  an  intermediate 
component  of  reflection  observed  in  flat-panel  displays  or 
haze,  caused  by  local  light  scattering.  The  spread  of  the 
specular  peak  shows  the  difficulty  in  defining  specifically 
the  specular  and  haze  regions.  We  note  that  the  diffrise 
component  follows  a  quasi -Lambertian  profile  with  about 
the  same  fall-off  rate  with  off-axis  angle  found  for  the 
CRT  models.  We  speculate  that  the  several  discontinu¬ 
ities  are  associated  with  refractive  index  changes  of  the 
different  layers  of  the  OLED  stack.  From  the  simulation 
results,  it  is  possible  to  compute  a  haze  ratio  associated 
with  the  full  width  at  half  maximum  of  about  2°,  How¬ 
ever,  as  described  in  experimental  studies,  the  value  of  the 
haze  ratio  is  very  sensitive  to  the  cut-off  limit.  In  Figs.  5 
and  6,  we  show  a  comparison  of  the  simulated  reflected 
beam  spot  profiles  for  the  CRT  and  OLED  model.  The 
significant  diffuse  reflections  off  the  CRT  phosphor  layer 
cause  blur.  For  the  OLED  model,  the  mirror-like  appear¬ 
ance  of  the  reflected  spot  is  consistent  with  a  small  diffuse 
component. 


IV.  Discussion 

We  present  a  novel  method  based  on  Monte  Carlo  tech¬ 
niques  to  calculate  the  bidirectional  reflectance  distribu¬ 
tion  function  of  emissive  displays  .  We  demonstrate  the 
validity  of  the  results  using  simple  models  for  which  the 
exact  analytical  solution  is  known,  and  by  comparing  the 
results  with  experimental  measurements.  We  conclude 
that  for  both  emissive  display  technologies,  the  reflections 
are  slightly  more  forward-peaked  than  the  predictions  of 
Lambertian  profiles.  We  show  that  the  BRDF  of  OLED 
displays  based  on  polymer  films  backed  with  a  metallic 
reflective  electrode  is  highly  specular.  We  report  for  the 
first  time  the  complete  BRDF  signature  of  OLEDs,  show¬ 
ing  that  the  haze  component  is  small,  with  a  characteristic 
width  of  2®.  This  intermediate  component  of  reflection 
complicates  the  measurement  of  the  specular  peak  since 
luminance  meters  having  slightly  different  apertures  will 
produce  a  specular  reflection  component  with  large  er¬ 
rors.  This  modelling  tool  will  be  used  to  demonstrate  im¬ 
provements  of  novel  designs  for  achieving  high-contrast 
in  high-illuminance  environments  without  the  need  to  fab¬ 
ricate  test  devices,  or  use  expensive  instrumentation.  We 
are  investigating  the  effect  of  the  optical  properties  and 
absorptive  layers  on  the  reflectance  properties  of  OLEDs, 
and  the  BRDF  signatures  of  liquid  crystal  displays. 

References 

[11  A.  Badano.  Image  Quality  Degradation  by  Light  Scattering  Pro¬ 
cesses  in  High  Performance  Display  Devices  for  Medical  Imag¬ 
ing.  PhD  thesis,  University  of  Michigan,  1999. 

[2]  A.  Badano.  Monte  Carlo  modeling  of  the  bidirectional  re¬ 
flectance  of  display  devices.  In  preparation^  2001 . 

[3]  A.  Badano,  M.  J.  Flynn,  E.  Muka,  K.  Compton,  and  T.  Mon- 
sees.  The  veiling  glare  point-spread  function  of  medical  imaging 
monitors.  In  Proceedings  of  the  SPIE  Medical  Imaging  Display 
Conference^  1999. 

[4]  A.  Badano  and  J,  Kanicki.  Monte  carlo  analysis  of  the  spec¬ 
tral  photon  emission  and  extraction  efficiency  of  organic  light- 
emitting  devices.  Journal  of  Applied  Physics,  90(4),  2001 . 

[5]  M.  E.  Becker.  Evaluation  and  characterization  of  display  re¬ 
flectance.  Displays,  19:35-54,  1998. 

[6]  H.  Blume.  In  Proceedings  of  the  SPIE  Medical  Imaging  Confer¬ 
ence,  volume  4323-07, 2001. 

[7]  M.  Bom  and  E.  Wolf.  Principles  of  Optics.  3rd  revised  edition, 
1965. 

[8]  J.  Delacour,  S.  Ungar,  G.  Mathieu,  G.  Hasna,  P.  Martinez,  and 
J.-C.  Roche.  Front  panel  engineering  with  CAD  simulation  tool. 
In  Proceedings  of  the  SPIE  Flat  Panel  Display  Technology  and 
Display  Metrology  Conference,  1999. 

[9]  M.  Elias,  L.  Simonot,  and  M.  Menu.  Bidirectional  reflectance  of 
a  diffuse  background  covered  by  a  partly  absorbing  layer.  Optics 
Communications,  191:1-7,  200 1 . 

[10]  Y.  He,  S.  Gong,  R.  Hattori,  and  J.  Kanicki.  Applied  Physics  Let¬ 
ters,  74, 1999. 

[11]  E.  F.  Kelley.  Display  reflectance  model  based  on  BRDF.  Dis¬ 
plays,  19:27-34, 1998. 


MONTE  CARLO  SIMULATION  OF  SPECTRAL  PHOTON  EMISSION  OF  THE  ORGANIC  POLYMER 

LIGHT-EMITTING  DEVICES 

Shu-ien  Lee,  *  Aldo  Badano,*^^  Yongtaek  Hong,®^  and  Jerzy  Kanicki®^ 

Organic  &  Molecular  Electronics  Laboratory,  Department  of  Electrical  Engineering  and  Computer  Science,  The 
University  of  Michigan  -Ann  Arbor,  Ml  48109 
Macromolecular  Science  and  Engineering  Center,  The  University  of  Michigan  —  Ann  Arbor 
U  S.  Food  and  Drug  Administration,  Center  for  Devices  and  Radiological  Health,  Rockville,  MD. 

Abstract 

We  reported  a  Monte  Carlo  method  for  modeling  the  light  transport  phenomena  in  the  organic  polymer 
light-emitting  devices  (OP-LEDs)  fabricated  on  plastic  substrates.  The  OP-LED  consists  of  a  combination  of  a  hole 
transport  layer  (HTL)  of  PEDOT/PSS  and  a  light-emitting  layer  of  Dow  Chemical  green  B  or  red  B  organic 
polymer,  Fig.  1.  In  this  simulation  we  assumed  a  point  light  source  having  photon  emission  spectra  represented  by 
photoluminescence  (PL)  spectra  of  the  organic  polymers.  This  simulation  method  describes  the  fate  of  photons 
through  multiple  scattering  events  determined  by  the  wavelength-dependent  material  optical  properties  in  a  3-D 
Cartesian  geometry,  thus  considering  the  effects  of  refraction  at  different  interfaces,  back-reflection,  and  absorption 
within  the  polymer  layers.  The  absorption  coefficients,  used  in  this  simulation  work  were  measured  by  a 
combination  of  the  UV-visible  absorption  and  the  photothermal  deflection  spectroscopy,  and  the  refractive  indices, 
n(\),  were  obtained  by  spectroscopic  ellipsometry. 

We  apply  this  method  to  analyze  the  light  wavelength  distribution  and  extraction  efficiency.  We  found  that 
the  simulated  light  emission  spectra  of  the  green  and  red  light-emitting  devices  are  very  similar  to  the  measured  PL 
spectra.  Fig.  2,  We  also  established  that  the  calculated  extraction  efficiency  for  the  red  (T|ext=21.71%)  and  green  (r|ext 
=21.70%)  OP-LEDs  are  approximately  the  same.  We  further  investigated  the  light  emission  angular  distribution  of 
the  OP-LEDs,  and  found  that  the  angular  distribution  is  rather  a  pseudo-Lambertian  than  Lambertian  with  increased 
intensity  in  the  forward  direction,  Fig.  3.  It  was  found  that  the  light  refraction  at  different  interfaces  leads  to  this 
pseudo-Lambertian  angular  distribution,  while  the  absorption  and  the  back-reflection  have  minor  effects.  Finally,  it 
is  clear  from  our  simulation  that  the  waveguiding  effect  cannot  be  neglected  in  the  OP-LEDs  and  its  contribution 
will  be  dependent  on  substrate  thickness  and  refractive  index  values. 
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Fig.  1:  The  cross-section  of  the  device  structure  used 
in  this  simulation  work  fnot  to  scaleL 
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Fig.  2:  The  simulated  light  output  spectra  of  the  OP- 
LEDs.  The  PL  spectra  for  green  B  and  red  B  organic 
polymers  are  also  shown  in  this  figure. 


Fig.  3:  The  simulated  light  emission  angular 
distribution  of  the  OP-LEDs. 
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CONCLUSIONS 

A  Monte  Carlo  method  was  uesd  for  modeling  the  light  transport  phenomena  of  the  organic 
polymer  light  emitting  devices 

The  caculated  extraction  efficiencies  for  the  red  and  green  polymer  devices  are  ^  21.7% 
The  simulated  output  angular  distribution  is  a  quasi  Lambertian  angular  distribution  with 
forward  peeked  intensity 

Further  research  is  needed  to  qualify  this  simulation  method  for  optimization  of  the  OP-LED 
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We  report  on  a  Monte  Carlo  method  for  modeling  light  transport  phenomena  in  multilayer  organic 
polymer  light-emitting  devices  on  plastic  flexible  substrates.  The  method  allows  modeling  of 
Cartesian  geometrical  structures  describing  the  fate  of  photons  through  multiple  scattering  events 
determined  by  the  wavelength-dependent  material  optical  properties.  We  apply  the  method  to 
analyze  the  wavelength  distribution  of  emitted  light  spectra.  We  find  that  for  all  organic  polymers 
considered,  the  light  emission  is  slightly  shifted  toward  the  longer  wavelengths,  and  that  this  shift 
is  maximum  for  light  emissions  with  peaks  around  530  nm.  The  photon  extraction  efficiency  is 
higher  (0.430)  for  organic  polymers  emitting  in  the  longer  wavelengths,  while  the  photon  absorbed 
fraction  is  higher  (0.676)  for  spectra  with  a  maximum  in  the  short  wavelengths.  ©  2001  American 
Institute  of  Physics,  [DOI:  10.1063/1.1385571] 


I.  INTRODUCTION 

The  increase  in  the  luminous  efficiency  of  organic  light- 
emitting  devices  (OLEDs)  remains  the  focus  of  many  efforts 
centered  mostly  on  controlling  charge  transport  and  increas¬ 
ing  the  photoluminescent  efficiency  of  organic  materials. 
Currently,  the  analysis  of  the  optical  transport  processes  in 
OLEDs  has  obtained  much  attention.  It  is  known  that  only 
about  one  fifth  of  the  light  generated  in  an  OLED  is  emitted 
through  the  top  surface  of  a  simple  OLED  structure  on  a 
glass  substrate.^  The  optical  losses  are  related  to  absorption 
in  the  organic  material  and  edge  emission  of  waveguided 
modes.  For  a  typical  organic  semiconductor  with  refractive 
index  of  2.0,  the  escape  cone  angle  into  a  substrate  with  an 
index  of  1.5  is  48.6°  which  corresponds  to  a  solid  angle  of 
only  ^^6%.  Several  attempts  have  been  made  to  reduce  the 
waveguided  light  fraction  by  using  structured  surfaces^""^  and 
films, ^  and  spherical  scatterers  within  the  organic  film.^  In 
addition  to  affecting  the  efficiency  of  the  OLED,  photon 
transport  processes  influence  the  spectral  and  angular  distri¬ 
butions  of  the  emitted  fight  intensity.  Both  distributions  are 
important  design  parameters  for  achieving  full-color  efficient 
OLEDs  for  flat-panel  display  applications. 

After  the  simple  analysis  of  OLED  efficiency  presented 
in  Ref.  1,  several  groups  have  investigated  the  effect  of  fight 
transport  in  multilayer  structures. The  method  we  present 
in  this  article  is  based  on  a  Monte  Carlo  (MC)  approach. 

II.  SIMULATION  METHOD 

The  MC  method  makes  use  of  the  generation  of  photons 
with  random  direction  according  to  a  distribution  function 
describing  the  nature  of  the  fight  emission.  In  this  analysis, 

“^Current  address:  U.S.  Food  and  Drug  Administration,  Center  for  Devices 
and  Radiological  Health,  12720  Twinbrook  Parkway,  Rockville,  MD 
20857;  Electronic  mail:  agb@cdrh.fda.gov 
'’taectronic  mail:  kanicki@eecs.umich.edu 
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the  fight  source  within  the  organic  polymer  layer  is  consid¬ 
ered  isotropic  from  a  single  point  situated  in  the  center  of  the 
device  (see  Fig.  1).  To  obtain  an  isotropic  distribution  of  the 
directional  cosines,  we  sample  the  three  directional  cosines 
that  define  the  photon  direction  according  to 

VT-^sin  277^2 » 

Cy  =  Vl“^lCOS  2  77^2  . 

where  and  ^2  uniformly  sampled  in  [0,1).  The  energy 
of  the  photon  source  is  defined  by  a  table  that  corresponds  to 
a  specific  spectral  fight  emission.  The  initial  photon  polariza¬ 
tion  vector  is  sampled  uniformly  in  the  477  space,  therefore 
assuming  unpolaiized  fight  emission.  The  photon  histories 
are  then  followed  through  a  sequence  of  interactions  that 
includes  absorption  and  Fresnel  refraction.  A  unique  advan¬ 
tage  of  this  simulation  method  is  its  ability  to  model  bulk 
absorption  events,  thin-film  coatings  and  rough  surfaces, 
while  keeping  track  of  the  photon  polarization  state.  Bulk 
absorption  is  determined  by  sampling  the  probability  of  a 
photon  being  absorbed  after  a  path  of  length  /  by  the  expo¬ 
nential  law 

where  ptabi^)  is  the  wavelength-dependent  linear  absorption 
coefficient.  At  the  optical  boundaries,  an  analysis  is  per¬ 
formed  depending  on  the  surface  type  and  material  proper¬ 
ties  using  Fresnel’s  equations  and  considering  the  polariza¬ 
tion  of  the  incoming  photon.  When  the  film  thickness  is 
comparable  to  the  photon  wavelength,  we  use  modified 
Fresnel  coefficients  to  describe  the  interference  effects  of 
optically  thin  films.  The  reflection  and  transmission  coeffi¬ 
cients  are  then  interpreted  as  probabilities.  The  simulation 
outcome  is  calculated  by  a  statistical  average  of  the  fate  of 
all  histories  according  to  the  desired  quantity  to  be  evaluated 
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FIG.  1.  Schematic  cross  section  of  the  OLED  structure  used  for  the  Monte 
Carlo  simulations,  showing  the  location  of  the  centered  point  light  source  at 
the  interface  between  the  hole-transporting  and  emissive  organic  polymer 
layers. 


for  each  experiment.  Possible  reporting  options  include  the 
angular  and  spectral  distributions  of  the  emitted  photons,  the 
point-spread  function,  the  specular  and  diffuse  reflection  co- 
efflcients,  and  a  summary  of  scattering  event  statistics. 

From  the  light  emission  at  the  luminescent  center  until 
the  photons  emerge,  multiple  scattering  events  take  place 
within  the  multilayer  OLED  structure.  For  the  purpose  of  our 
analysis,  we  define  the  device  external  quantum  efficiency  as 

Vin  Vpe  ’ 

where  is  the  intrinsic  quantum  efficiency  related  to  car¬ 
rier  recombination  and  photoluminescent  fraction,  and  77^^^  is 
the  photon  extraction  efficiency.  We  introduce  77^^  to  repre¬ 
sent  the  probability  that  a  photon  generated  at  the  lumines¬ 
cent  center  within  the  OLED,  emerges  through  the  front  sur¬ 
face  of  the  device  (through  the  transparent  electrode), 
thereby  contributing  to  luminance.  The  77^,^  depends  strongly 
on  the  device  structure  and  on  the  material  and  surface  prop¬ 
erties,  and  is  always  less  than  unity  due  to  light  absorption, 
waveguiding,  and  edge  emissions.  We  can  summarize  the 
relevant  physical  processes  that  occur  as 

Vpe  ^  Vwa  Vab  Vtr'> 

where  77^^  is  the  fraction  of  photons  that  are  waveguided 
within  the  structure  and  exit  through  the  device  edges,  77^/,  is 
the  absorbed  fraction,  and  77^^  is  the  fraction  transmitted 
through  the  metallic  electrode  deposited  in  the  side  opposite 
the  direction  of  the  desired  light  emission.  The  top  metallic 
cathode  electrode  of  all  the  OLED  structures  modeled  in  this 
article  is  an  aluminum  layer  deposited  by  vacuum 
evaporation.^^  We  considered  77,^= 0  in  all  simulations  pre¬ 
sented  in  this  article.  In  this  work,  we  neglected  the  photo¬ 
luminescence  quenching  due  to  polymer  composition  or 
blend  variations  and  the  presence  of  carrier  flow  within  the 
OLED.  Electric  field  induced  photoluminescence  (PL) 
quenching  in  conjugated  polymers,  which  is  caused  by  exci- 
ton  dissociation,  is  well  known, but  cannot  be  implemented 
easily  in  this  calculation.  It  should  be  noted  that  this  PL 
quenching  is  not  important  at  the  OLED  operating  point. 

111.  RESULTS  AND  DISCUSSION 

We  measured  the  index  of  refraction,  absorption,  PL,  and 
electro-luminescence  (EL)  of  three  organic  polymers  (A,  B, 
and  C)  with  peak  emissions  in  a  different  region  of  the  vis¬ 
ible  spectrum.  The  refractive  index  (considered  equal  for  the 
three  organic  polymers)  and  absorption  characteristics  are 
shown  in  Fig.  2.  We  used  the  PL  spectrum  as  the  photon 


FIG.  2.  Refractive  index  and  normalized  absorption  coefficient  of  organic 
polymers  A,  B,  and  C  modeled  in  this  work. 


energy  distribution  at  the  source  for  the  MC  histories.  The 
OLED  structure  modeled  in  this  article  (see  Fig.  1)  is  a  het¬ 
erostructure  OLED  described  by  He  et  alP  as  one  with  an 
aluminum  cathode  electrode  and  a  transparent  anode  ITO 
electrode  (160  nm,  refractive  index  of  1.8).  The  organic 
polymer  film  thickness  used  was  200  nm.  The  transparent 
substrate  index  of  refraction  was  1.5  and  its  thickness  was 
900  fjLm. 

In  Fig.  3,  we  present  results  of  the  simulated  light  (SL) 
emission  wavelength  distribution.  For  the  three  PL  spectra 
considered,  the  measured  EL  spectrum  is  only  slightly 
shifted  toward  the  longer  wavelengths.  Our  MC  simulation 
results  are  consistent  with  this  trend.  By  computing  sepa¬ 
rately  the  simulated  emission  from  a  device  structure  having 
absorption  in  the  organic  polymer  film,  and  from  another 
with  no  absorption  but  with  a  thin-film  transparent  layer  be¬ 
tween  the  organic  material  and  the  substrate,  we  have  con¬ 
firmed  that  the  decrease  in  power  efficiency  of  the  shorter 
wavelength  range  is  associated  with  absorption  in  the  organic 
material,  while  the  increase  in  strength  at  longer  wavelengths 
is  caused  by  interference  effects  associated  with  the  transpar¬ 
ent  conductive  coating.  Our  simulation  method  can  correctly 
predict  the  shift  in  maximum  wavelength  of  each  EL  spec¬ 
trum  with  respect  to  the  PL  spectrum  of  each  material.  It 
should  be  noticed  that  the  maximum  of  the  SL  spectrum  of 
each  material  is  located  very  close  to  the  corresponding 
maximum  of  the  measured  EL  spectrum.  However,  espe¬ 
cially  in  the  longer  wavelengths,  a  discrepancy  exists  be¬ 
tween  the  EL  and  the  MC  calculated  SL  that  cannot  be  ex¬ 
plained  in  terms  of  optical  transport  phenomena.  The 
discrepancy  between  the  SL  and  EL  suggests  that  the  spectral 
light  emission  at  the  source  might  not  be  completely  accu¬ 
rate.  For  our  MC  simulations,  we  have  assumed  that  the 
spectrum  at  the  source  is  equal  to  the  PL  spectrum  measured 
for  each  polymer.  The  actual  light  emission  in  OLEDs  could 
differ  from  that  of  the  PL,  reflecting  the  geometrical  and 
structural  differences  in  both  the  solid-state  films  in  the 
OLED  and  in  the  device  structure.  Photoluminescence 
quenching  due  to  the  presence  of  carriers  might  also  affect 
the  spectral  emission  by  modifying  the  ratio  between  radia¬ 
tive  and  nonradiative  recombination. 
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FIG.  3.  Photoluminescent,  electroluminescent,  and  simulated  light  emission 
spectra  of  organic  polymers  A,  B,  and  C. 


In  general,  it  is  expected,  and  was  experimentally  con¬ 
firmed  for  our  OLEDs,  that  the  intensity  of  the  PL  spectrum 
decreases  with  increasing  applied  voltage  (e,g.,  with  increas¬ 
ing  current  density).  However,  this  decrease  is  not  very  large. 
Furthermore,  it  is  expected  that  the  light  emission  from  ex- 
cimers  can  contribute  to  the  EL  spectrum  at  longer  wave¬ 
lengths.  This  contribution  will  depend  on  the  polymer  struc¬ 
ture  and  will  be  significant  in  copolymers.  We  have  checked 
experimentally  that  the  PL  spectra  shape  is  independent  of 
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FIG.  4.  Extraction  efficiency  rjpg ,  waveguiding  7jy,a ,  and  absorption 
for  different  transparent  conductive  oxide  electrode  models:  (a)  effect  of 
refractive  index,  and  (b)  effect  of  thin-film  thickness.  The  vertical  line  rep¬ 
resents  the  default  values  for  the  refractive  index  and  thickness  of  the  trans¬ 
parent  electrode  used  for  the  calculations  of  the  light  emission  spectra 
shown  in  Fig.  3.  Curves  are  shown  for  each  of  the  three  polymers  modeled 
in  this  article  (A,  B,  and  C). 

film  thickness.  We  have  also  verified  through  simulation  that 
the  film  thickness  of  the  EL  layer  and  the  location  of  the 
light-emitting  source  within  the  OLED  structure  do  not  affect 
either  the  EL  spectra  shape  or  the  peak  maximum  location. 

The  extraction  efficiency  is  also  affected  by  the 
wavelength  distribution  of  the  photon  source.  For  the  case  of 
a  transparent  electrode  having  a  thickness  of  160  nm  and  a 
refractive  index  of  1.80,  97^^  is  0.314  for  polymer  A,  0.334 
for  polymer  B,  and  0.430  for  polymer  C.  The  absorbed  frac¬ 
tion  rjab  is  0.676,  0.655,  and  0.553,  and  the  wave-guided 
fraction  77^^  is  0.001,  0.011,  and  0.017,  respectively,  for 
polymers  A,  B,  and  C.  The  low  is  caused  by  high  ab¬ 
sorption  in  the  organic  film.  In  Fig.  4,  we  present  results 
showing  the  change  in  the  values  of  77^^ ,  77^^ ,  and  77^^  with 
the  index  of  refraction  (thickness)  of  the  transparent  conduct¬ 
ing  oxide  (TCO)  electrode  covering  the  range  from  1.1  to  2.7 
(10-200  nm).  We  find  that  both  the  thickness  and  the  refrac¬ 
tive  index  of  the  transparent  electrode  have  a  minor  effect  on 
the  shift  of  the  SL  spectrum  with  respect  to  the  PL  emission. 
We  obtained  a  maximum  77^^  for  a  TCO  thickness  and  a 
refractive  index  of  about  170  nm  and  1.85,  respectively. 
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FIG.  5.  Extraction  efficiency  77^^ ,  waveguiding  7^.^ ,  and  absorption  7^^  as 
a  function  of  OLED  size  for  polymer  B.  Similar  curves  were  obtained  for 
polymers  A  and  C.  The  vertical  line  represents  the  default  device  size  used 
throughout  this  article. 

The  waveguide  modes  are  determined  by  the  geometry 
of  the  OLED  stack  that  defines  a  total  internal  reflection 
(TIR)  angle  ( a™)  with  respect  to  the  device  plane,  beyond 
which  all  photons  are  emitted  through  the  edges.  Since  our 
source  of  photons  is  isotropic,  we  expect  lower  for 
laiger  devices  due  to  (1)  an  increased  probability  of  reflec¬ 
tion  and  scattering  going  into  the  solid  angle  defined  by 
or^iR,  and  (2)  increased  absorption  in  the  organic  polymer 
layer.  The  simulation  result  presented  in  Fig.  5  confirms  this 
assumption.  For  OLED  sizes  larger  than  50  /zm,  77^^  steadily 
increases  because  more  photons  impinge  directly  (as  a  first 
interaction)  into  the  substrate.  For  small  device  sizes  below  2 
nun,  the  absorbed  fraction  77^^  remains  constant  while  the 
waveguiding  decreases  as  the  photon  extraction  efficiency 
increases.  When  the  device  size  decreases  below  100  /im, 
lateral  boundary  conditions  become  important  and  edge  ef¬ 
fects  are  observed.  In  our  simulations,  no  additional  structure 
is  present  at  the  edges  of  the  device  (i.e.,  the  sides  are  con¬ 
sidered  to  be  smooth  surfaces  in  contact  with  air).  Photons 
reflected  from  the  sides  of  the  substrate  at  larger  angles  with 
respect  to  the  surface  normal  have  a  larger  probability  of 
exiting  the  device  through  the  top  surface,  thereby  contribut¬ 
ing  to  the  extraction  efficiency  and  to  the  luminance.  For 
large  device  sizes,  77^^  and  77^^  increase  while  the  waveguid¬ 
ing  continues  to  decrease  at  a  faster  rate.  The  results  also 
confirm  that  since  most  of  the  photon  waveguiding  occurs  in 
the  transparent  substrate,  the  absorption  in  the  organic  poly¬ 
mer  film  has  a  minor  effect  on  reducing  the  light  edge  emis¬ 
sion  and  increasing  the  photon  extraction  efficiency. 


iV.  CONCLUSION 

In  this  article,  we  have  presented  a  method  to  model 
light  transport  processes  in  OLEDs  based  on  Monte  Carlo 
techniques.  For  all  organic  polymers  considered  in  this  woiic, 
our  results  show  that  the  simulated  light  emission  is  shifted 
toward  the  longer  wavelengths,  consistent  with  experimental 
measurements  of  EL.  We  also  showed  that  the  photon  extrac¬ 
tion  efficiency  is  reduced  by  light  absorption  and  waveguid¬ 
ing.  Reduction  of  the  OLED  waveguided  fraction  can  be 
obtained  by  using  structured  substrates  or  by  tuning  the  re¬ 
fractive  index  of  the  layers  to  maximize  the  transmission  and 
reduce  the  total  internal  reflection  at  each  interface.  Of  par¬ 
ticular  importance  to  this  approach  is  the  interface  between 
the  organic  polymers  and  the  TCO/substrate,  as  well  as  de¬ 
tails  of  the  structure  and  index  of  refraction  of  the  layers  in 
plastic  substrates.  The  presence  of  optimized  thin-film  coat¬ 
ings  in  a  plastic  substrate  could  contribute  to  a  decrease  in 
internally  reflected  trapped  photons.  A  combination  of  these 
two  approaches  could  result  in  optimal  OLED  structures 
with  high  photon  extraction  efficiency. 
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